Using Expression and Genotype to Predict Drug Response in Yeast

Personalized, or genomic, medicine entails tailoring pharmacological therapies according to individual genetic variation at genomic loci encoding proteins in drug-response pathways. It has been previously shown that steady-state mRNA expression can be used to predict the drug response (i.e., sensitivity or resistance) of non-genotyped mammalian cancer cell lines to chemotherapeutic agents. In a real-world setting, clinicians would have access to both steady-state expression levels of patient tissue(s) and a patient's genotypic profile, and yet the predictive power of transcripts versus markers is not well understood. We have previously shown that a collection of genotyped and expression-profiled yeast strains can provide a model for personalized medicine. Here we compare the predictive power of 6,229 steady-state mRNA transcript levels and 2,894 genotyped markers using a pattern recognition algorithm. We were able to predict with over 70% accuracy the drug sensitivity of 104 individual genotyped yeast strains derived from a cross between a laboratory strain and a wild isolate. We observe that, independently of drug mechanism of action, both transcripts and markers can accurately predict drug response. Marker-based prediction is usually more accurate than transcript-based prediction, likely reflecting the genetic determination of gene expression in this cross

Amino Acid Metabolic Origin as an Evolutionary Influence on Protein Sequence in Yeast

The metabolic cycle of Saccharomyces cerevisiae consists of alternating oxidative (respiration) and reductive (glycolysis) energy-yielding reactions. The intracellular concentrations of amino acid precursors generated by these reactions oscillate accordingly, attaining maximal concentration during the middle of their respective yeast metabolic cycle phases. Typically, the amino acids themselves are most abundant at the end of their precursor’s phase. We show that this metabolic cycling has likely biased the amino acid composition of proteins across the S. cerevisiae genome. In particular, we observed that the metabolic source of amino acids is the single most important source of variation in the amino acid compositions of functionally related proteins and that this signal appears only in (facultative) organisms using both oxidative and reductive metabolism. Periodically expressed proteins are enriched for amino acids generated in the preceding phase of the metabolic cycle. Proteins expressed during the oxidative phase contain more glycolysis-derived amino acids, whereas proteins expressed during the reductive phase contain more respiration-derived amino acids. Rare amino acids (e.g., tryptophan) are greatly overrepresented or underrepresented, relative to the proteomic average, in periodically expressed proteins, whereas common amino acids vary by a few percent. Genome-wide, we infer that 20,000 to 60,000 residues have been modified by this previously unappreciated pressure. This trend is strongest in ancient proteins, suggesting that oscillating endogenous amino acid availability exerted genome-wide selective pressure on protein sequences across evolutionary time

Genetic Architecture of Highly Complex Chemical Resistance Traits across Four Yeast Strains

Many questions about the genetic basis of complex traits remain unanswered. This is in part due to the low statistical power of traditional genetic mapping studies. We used a statistically powerful approach, extreme QTL mapping (X-QTL), to identify the genetic basis of resistance to 13 chemicals in all 6 pairwise crosses of four ecologically and genetically diverse yeast strains, and we detected a total of more than 800 loci. We found that the number of loci detected in each experiment was primarily a function of the trait (explaining 46% of the variance) rather than the cross (11%), suggesting that the level of genetic complexity is a consistent property of a trait across different genetic backgrounds. Further, we observed that most loci had trait-specific effects, although a small number of loci with effects in many conditions were identified. We used the patterns of resistance and susceptibility alleles in the four parent strains to make inferences about the allele frequency spectrum of functional variants. We also observed evidence of more complex allelic series at a number of loci, as well as strain-specific signatures of selection. These results improve our understanding of complex traits in yeast and have implications for study design in other organisms

Genetic Basis of Hidden Phenotypic Variation Revealed by Increased Translational Readthrough in Yeast

Eukaryotic release factors 1 and 3, encoded by SUP45 and SUP35, respectively, in Saccharomyces cerevisiae, are required for translation termination. Recent studies have shown that, besides these two key factors, several genetic and epigenetic mechanisms modulate the efficiency of translation termination. These mechanisms, through modifying translation termination fidelity, were shown to affect various cellular processes, such as mRNA degradation, and in some cases could confer a beneficial phenotype to the cell. The most studied example of such a mechanism is [PSI+], the prion conformation of Sup35p, which can have pleiotropic effects on growth that vary among different yeast strains. However, genetic loci underlying such readthrough-dependent, background-specific phenotypes have yet to be identified. Here, we used sup35C653R, a partial loss-of-function allele of the SUP35 previously shown to increase readthrough of stop codons and recapitulate some [PSI+]-dependent phenotypes, to study the genetic basis of phenotypes revealed by increased translational readthrough in two divergent yeast strains: BY4724 (a laboratory strain) and RM11_1a (a wine strain). We first identified growth conditions in which increased readthrough of stop codons by sup35C653R resulted in different growth responses between these two strains. We then used a recently developed linkage mapping technique, extreme QTL mapping (X-QTL), to identify readthrough-dependent loci for the observed growth differences. We further showed that variation in SKY1, an SR protein kinase, underlies a readthrough-dependent locus observed for growth on diamide and hydrogen peroxide. We found that the allelic state of SKY1 interacts with readthrough level and the genetic background to determine growth rate in these two conditions

Variations in Stress Sensitivity and Genomic Expression in Diverse S. cerevisiae Isolates

Interactions between an organism and its environment can significantly influence phenotypic evolution. A first step toward understanding this process is to characterize phenotypic diversity within and between populations. We explored the phenotypic variation in stress sensitivity and genomic expression in a large panel of Saccharomyces strains collected from diverse environments. We measured the sensitivity of 52 strains to 14 environmental conditions, compared genomic expression in 18 strains, and identified gene copy-number variations in six of these isolates. Our results demonstrate a large degree of phenotypic variation in stress sensitivity and gene expression. Analysis of these datasets reveals relationships between strains from similar niches, suggests common and unique features of yeast habitats, and implicates genes whose variable expression is linked to stress resistance. Using a simple metric to suggest cases of selection, we found that strains collected from oak exudates are phenotypically more similar than expected based on their genetic diversity, while sake and vineyard isolates display more diverse phenotypes than expected under a neutral model. We also show that the laboratory strain S288c is phenotypically distinct from all of the other strains studied here, in terms of stress sensitivity, gene expression, Ty copy number, mitochondrial content, and gene-dosage control. These results highlight the value of understanding the genetic basis of phenotypic variation and raise caution about using laboratory strains for comparative genomics

Using Stochastic Causal Trees to Augment Bayesian Networks for Modeling eQTL Datasets

<p>Abstract</p> <p>Background</p> <p>The combination of genotypic and genome-wide expression data arising from segregating populations offers an unprecedented opportunity to model and dissect complex phenotypes. The immense potential offered by these data derives from the fact that genotypic variation is the sole source of perturbation and can therefore be used to reconcile changes in gene expression programs with the parental genotypes. To date, several methodologies have been developed for modeling eQTL data. These methods generally leverage genotypic data to resolve causal relationships among gene pairs implicated as associates in the expression data. In particular, leading studies have augmented Bayesian networks with genotypic data, providing a powerful framework for learning and modeling causal relationships. While these initial efforts have provided promising results, one major drawback associated with these methods is that they are generally limited to resolving causal orderings for transcripts most proximal to the genomic loci. In this manuscript, we present a probabilistic method capable of learning the causal relationships between transcripts at all levels in the network. We use the information provided by our method as a prior for Bayesian network structure learning, resulting in enhanced performance for gene network reconstruction.</p> <p>Results</p> <p>Using established protocols to synthesize eQTL networks and corresponding data, we show that our method achieves improved performance over existing leading methods. For the goal of gene network reconstruction, our method achieves improvements in recall ranging from 20% to 90% across a broad range of precision levels and for datasets of varying sample sizes. Additionally, we show that the learned networks can be utilized for expression quantitative trait loci mapping, resulting in upwards of 10-fold increases in recall over traditional univariate mapping.</p> <p>Conclusions</p> <p>Using the information from our method as a prior for Bayesian network structure learning yields large improvements in accuracy for the tasks of gene network reconstruction and expression quantitative trait loci mapping. In particular, our method is effective for establishing causal relationships between transcripts located both proximally and distally from genomic loci.</p

Joint Genetic Analysis of Gene Expression Data with Inferred Cellular Phenotypes

Even within a defined cell type, the expression level of a gene differs in individual samples. The effects of genotype, measured factors such as environmental conditions, and their interactions have been explored in recent studies. Methods have also been developed to identify unmeasured intermediate factors that coherently influence transcript levels of multiple genes. Here, we show how to bring these two approaches together and analyse genetic effects in the context of inferred determinants of gene expression. We use a sparse factor analysis model to infer hidden factors, which we treat as intermediate cellular phenotypes that in turn affect gene expression in a yeast dataset. We find that the inferred phenotypes are associated with locus genotypes and environmental conditions and can explain genetic associations to genes in trans. For the first time, we consider and find interactions between genotype and intermediate phenotypes inferred from gene expression levels, complementing and extending established results

Methodological approaches in application of synthetic lethality screening towards anticancer therapy

A promising direction in the development of selective less toxic cancer drugs is the usage of synthetic lethality concept. The availability of large-scale synthetic low-molecular-weight chemical libraries has allowed HTS for compounds synergistic lethal with defined human cancer aberrations in activated oncogenes or tumour suppressor genes. The search for synthetic lethal chemicals in human/mouse tumour cells is greatly aided by a prior knowledge of relevant signalling and DNA repair pathways, allowing for educated guesses on the preferred potential therapeutic targets. The recent generation of human/rodents genome-wide siRNAs, and shRNA-expressing libraries, should further advance this more focused approach to cancer drug discovery

Phenotypic Landscape of Saccharomyces cerevisiae during Wine Fermentation: Evidence for Origin-Dependent Metabolic Traits

The species Saccharomyces cerevisiae includes natural strains, clinical isolates, and a large number of strains used in human activities. The aim of this work was to investigate how the adaptation to a broad range of ecological niches may have selectively shaped the yeast metabolic network to generate specific phenotypes. Using 72 S. cerevisiae strains collected from various sources, we provide, for the first time, a population-scale picture of the fermentative metabolic traits found in the S. cerevisiae species under wine making conditions. Considerable phenotypic variation was found suggesting that this yeast employs diverse metabolic strategies to face environmental constraints. Several groups of strains can be distinguished from the entire population on the basis of specific traits. Strains accustomed to growing in the presence of high sugar concentrations, such as wine yeasts and strains obtained from fruits, were able to achieve fermentation, whereas natural yeasts isolated from “poor-sugar” environments, such as oak trees or plants, were not. Commercial wine yeasts clearly appeared as a subset of vineyard isolates, and were mainly differentiated by their fermentative performances as well as their low acetate production. Overall, the emergence of the origin-dependent properties of the strains provides evidence for a phenotypic evolution driven by environmental constraints and/or human selection within S. cerevisiae